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ABSTRACT 

Aims. To reconstruct the spectral shape of the quasar ionizing radiation in the extreme-UV range (1 Ryd < i5 < 10 
Ryd) from the analysis of narrow absorption lines (NAL) of the associated systems. 

Methods. Computational technique for inverse spectroscopic problems - Monte Carlo Inversion augmented by procedure 
of the spectral shape recovering and modified to account for the incomplete coverage of the light source. 
Results. The ionizing spectra responsible for the ionization structure of the NAL systems require an intensity depression 
at £ > 4 Ryd which is attributed to the Hell Lyman continuum opacity (r^°" ~1). A most likely source of this opacity 
is a quasar accretion disk wind. The corresponding column density of H I in the wind is estimated as a few times 10^® 
cm~^ . This amount of neutral hydrogen should cause a weak continuum depression at A < 912 A (rest-frame), and 
a broad and shallow absorption in H I Ly-a. If metallicity of the wind is high enough, other resonance lines of O vi, 
Nevi - Nevill, etc. are expected. In the analyzed QSO spectra we do observe broad (stretching over 1000s km s"^ ) and 
shallow (r ^ 1) absorption troughs of Hi Ly-Q and O vi AA1031, 1037 as well as continuum depressions at A < 912 A 
which correspond to A'^(Hl) ~5 x lO^'' cm~^ . Observational data available in both the UV and X-ray ranges suggest 
that at least ~50% of the quasar radiation passes through the gas opaque in the He ll Lyman continuum. This means 
that the outcoming ionizing spectrum should have a pronounced intensity break at _E > 4 Ryd with a depth of this 
break depending on the angle to the rotational axis of the accretion disk (the larger the angle the deeper the break). 
The QSO spectra with a discontinuity at 4 Ryd can influence the rate of the Hen reionization in the intergalactic 
medium and partly explain inhomogeneous (patchy) ionization structure of the intergalactic Hen observed at z ~ 3. 

Key words. Cosmology: observations - Line: formation - Line: profiles - Quasars: absorption lines 



1. Introduction in principle, impossible due to the intergalactic gas opacity 

in Hell Ly-a. 

5^ ; Spectral energy distribution (SED) of the quasar/ AGN out- Usually the gap in the SED between A = 912 A and soft 

. P_ , coming radiation in the extreme-UV range (EUV, A < 912 X-ray range is approximated by a simple power law, cx 

A) is important for the interpretation of both emission i/-". However, this part of the QSO continuum radiation is 

and absorption spectra since this range contains the ioniza- responsible for the ionization state of the so-called narrow 

tion thresholds of usually observed ions Cil-Civ, Nii-Nv, absorption Hne (NAL) systems - metal systems arising in 

O i-O VI, Si ii-Si IV. SED is accessible for direct measure- the gas located close to the quasar host galaxy - and can 

ments in the FUV-EUV ranges (rest-frame) for wavelengths be reconstructed from their analysis. 

A > 300 A (e.g., Zheng et al. 1997; Telfer et al. 2002; Scott The computational procedure aimed at restoring the 

et al. 2004) and at shorter wavelengths only in the X-ray spectral shape of the underlying ionizing radiation from 

range for A < 40 A (e.g., Piconcelli et al. 2005; Brocksopp et ions observed in optically thin metal absorption-line sys- 

al. 2006; Costantini et al. 2007). In the intermediate range tems is described in Agafonova et al. (2005, 2007). Applied 

(40 A < A < 300 A), direct measurements of the outcom- to the analysis of NAL systems identified in the spectrum 

ing radiation of QSOs with redshifts z < 2 are prevented of the quasar HE 0141-3932 with Zom= 1.8 (Reimers et 

by the Galactic absorption, and for QSOs with z > 2 are, al. 2005) this procedure resulted in a SED which could not 

be described with a single power law index a: the EUV 
Send offpr^nt requests to: S. A. Levshakov spectrum shows a step-like structure with a sharp intensity 
lev@astro.ioffe.rssi.ru break at A ~ 240 A (£, ~ 3.8 Ryd, shifted by ^12 A, or 
* Based on observations obtained at the VLT Kueyen tele- --16000 km s'^ redward from the Hell Ly-a edge 228 A, or 
scope (ESO, Paranal, Chile), and at the HST. 4 Ryd) and a much more slower intensity decay in the range 
** On leave from the loffe Physico- Technical Institute, St. 4 Ryd < i5 < 10 Ryd. Being attributed to the Lyman con- 
Petersburg, Russia tinuum absorption in He II, the depth of the break gives the 
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Hen column density 7V(Hell) = 0.7 x 10^* cm"^ . Similar 
spectra with breaks at ~4 Ryd corresponding to A'^(Hell) = 
(0.7 — 1.4) X 10^^ cm~^ were reconstructed from the inter- 
galactic absorption- line systems with Zabs ^1-8 (Fig. 20 
and Table 4 in Agafonova ct al. 2007). Since there are in- 
dications that a large part of the QSO/AGN emitted ra- 
diation remains at this redshift unprocessed by the IGM, 
we argued that these spectra may represent the energy dis- 
tribution of the outgoing quasar radiation. In the present 
paper, we continue the study of the intrinsic quasar SED 
using four associated (i.e. physically related to QSO) NAL 
systems identified in spectra of high-redshift (2 < < 3) 
QSOs. 

Concerning the selection of the NAL systems con- 
sidered in the paper the following is to be noted. The 
metal NAL systems are in general not rare - they are 
observed in 30%-50% of quasar spectra (Misawa et al. 
2007). However, for most of them the observed wavelength 
coverage and the ionization conditions in the absorbing 
gas result in the fact that at best only three highly ion- 
ized doublets CIVAA1548, 1551 A, NvAA1238, 1242 A, and 
O viAAlOSl, 1037 A arc available for the study. The column 
densities of these absorption lines can be reproduced with a 
wide variety of the ionizing spectra and, thus, do not allow 
to distinguish between specific SEDs. Another aggravating 
factor while working with NAL systems is that the esti- 
mation of the ion column densities is often hampered by 
effects like partial coverage of the background source of the 
continuum radiation and/or line blending. 

This explains the selection criteria which governed the 
search for the appropriate associated systems: (1) the pres- 
ence of many lines of different ions to obtain reasonable 
restrictions on the SED, (2) clear line profiles, and (3) the 
possibility to determine the covering factor with a suffi- 
ciently high accuracy. These requirements significantly re- 
stricted the number of suitable systems, but the loss in the 
quantity was compensated by the increased quality of the 
data which allowed us to clarify some additional aspects of 
the physical properties of the NAL systems. 

The paper is organized as follows. In Sect. 2 we re- 
construct SED of the intrinsic quasar EUV radiation from 
the analysis of quasar-related metal absorption-line systems 
identified in spectra of four QSOs with rcdshifts Zom = 2.2 
2.9. The results obtained are discussed in Sect. 3, and sum- 
marized in Sect. 4. 



2. Analysis of individual NAL systems 

Absorption systems arc analyzcxl by means of the Monte 
Carlo Inversion (MCI) procedure (Levshakov et al. 2000, 
2003; Agafonova et al. 2005, 2007) which is based on the 
assumption that all lines observed in the metal system are 
formed in the same gas with fluctuating density and ve- 
locity fields. The procedure's inputs are the observed line 
profiles and the ionization curves for each ion (hydrogen -|- 
metals) included in the analysis, the outputs - parameters 
of the gas density and velocity distributions and individ- 
ual metal abundances which are supposed to be constant 
within the absorber. The column densities of hydrogen and 
all metal ions are calculated as well. The ionization curves 
are computed with the photoionization code CLOUDY ver- 
sion 07.02.01 (last described by Fcrland ct al. 1998) which 
in turn uses as an input the adopted (trial) SED. 



For the present work the computational procedure was 
modified in a way to account for the incomplete coverage 
of the background light source. With C standing for the 
covering factor (fraction of the background source of the 
continuum radiation occulted by the absorbing cloud) and 
basing on principles of the geometrical optics one can easily 
determine the observed normalized intensity Ix within the 
absorption line profile as 

/A = Cexp(-TA) + (l-C). (1) 

The apparent optical depth t\ is calculated through Eq. (4) 
in Levshakov et al. (2000). The covering factor is assumed 
to be the same for the whole cloud, i.e. it does not depend 
on A (or, equivalently, on the radial velocity v). However, 
different covering factors arc allowed for different ions. 

It is well known that NAL systems with incomplete cov- 
erage often show A-depending covering factors which may 
be caused by the overlapping of clouds with different cov- 
ering factors or by the interaction of several physical pro- 
cesses. In these cases the estimation of accurate column 
densities becomes problematic and in the present study 
such systems are not used for the SED recovery. 

The SED of the radiation which ionizes the NAL sys- 
tem is reconstructed according to the iterative procedure 
described in Agafonova et al. (2005, 2007). The procedure 
includes (1) the parameterization of the spectral shape by 
means of a set of variables (called 'factors'), (2) the choice 
of a quantitative measure (called 'response') to evaluate 
goodness of a trial spectral shape, and (3) the estimation 
of a direction in the factor space which leads to a spectrum 
with the better goodness. Moving along this direction, we 
come from the initial spectral shape to the one with better 
characteristics concerning the fitting of the observed line 
intensities. Now the MCI calculations are carried out with 
this newly obtained ionizing spectrum and the whole pro- 
cedure is repeated till the optimal spectral shape is found, 
i.e. one which allows to reproduce the observed intensities 
of all lines without any physical inconsistencies. 

All calculations throughout the paper were performed 
with laboratory wavelength and oscillator strengths taken 
from Morton (2003) for A > 912 A and from Vcrner ct al. 
(1994) for A < 912 A. Solar abundances were taken from 
Asplund et al. (2004). Note that their solar abundances of 
neon, Ne/H = 7.5 x 10-^ and nitrogen, N/H = 6 x 10"^, 
are 1.4 times (0.15 dex) lower than those from Holweger 
(2001) which are used as default inputs in CLOUDY. 

2.1. NAL systems towards J 2233-606 

The quasar J 2233-606 (2:em= 2.238) was observed both 
at the VLT/UVES and HST/STIS with a wide wave- 
length coverage from 2300 A to ~ 10000 A (Cristiani & 
D'Odorico 2000). The spectrum reveals an associated com- 
plex consisting of several absorption systems extending over 
1000 km s~^ . Due to the wavelength coverage it becomes 
possible to identify not only the usually observed lines 
of the doublets C ivAA1548, 1551, N vAA1238, 1242, and 
OviAA1032, 1037, but also the lines of Ne viiiAA770, 780 
(albeit in a very noisy part of the STIS spectrum). 

The associated systems towards J 2233-606 have been 
studied twice - by Petitjean & Srianand (1999) and by 
Gabel et al. (2006). Both groups determined the column 
densities of different species, Ni, using velocity-dependent 
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Fig. 1. Hydrogen and metal absorption lines from the Zabs = 2.198 system towards J 2233-606 (solid-line histograms). The vertical 
axis is normalized intensities. The zero radial velocity is fixed at z = 2.1982. Synthetic profiles are plotted by the smooth curves. 
In panel Neviii, the dashed curve corresponds to the profile predicted for the power law ionizing spectrum oc v^^'^, and 
the smooth curve - to the profile predicted for the SED with a break at ~4 Ryd (Fig. 3). The observed profile of OivA787 A 
(dotted-line histogram) overplotted on O VI A1031 is blended with Ly-7 from the Zabs = 1.59 system in the range ~20-70 km s~^ . 



covering factors and after that compared iVj with the val- 
ues predicted by models where simple power law SEDs 
in the range from UV to X-ray regions were assumed. It 
was concluded that such models cannot account for the 
observed strong Neviii absorption and that an additional 
high-ionization absorber is needed to reproduce the Ne VIII 
line. In the present section we re-analyze two absorption 
systems from this complex using physical model of the as- 
sociated absorber as outlined above. Namely, gas in the 
absorbing cloud has varying density and velocity and the 
covering factor is assumed to be identical for all points in 
the line profile (no velocity dependence). However, the cov- 
ering factor can differ for individual ions since different ions 
trace different parts of the gas resulting in different solid 
angles subtended. 



2.1.1. System at 2:abs= 2.198 

Apart from metal doublets, this system exhibits also clear 
lines of hydrogen Ly-a and Ly-/? (Fig. 1). There are con- 
tinuum windows in the QSO spectrum at the expected po- 
sitions of C111A977 A and SiivA1393 A, the expected po- 
sition of O IVA787 A is blended with Ly-7 from the ^abs = 
1.59 system (O iv profile is overplotted on the panel with 
OviAlOSl A in Fig. 1). 

First trial calculations were performed with the power 
law SED, a = 1.5. The objective function included hydro- 
gen lines plus all lines of metal doublets mentioned above 
except for Neviii because of its too low S/N ratio. Since 
the system under consideration docs not show ionic lines of 
subsequent ionization stages, additional assumptions are re- 
quired to fix the mean ionization parameter U . For overso- 
lar metallicity which is reliably reproduced for this system 
independently on the model assumptions, the usual con- 
straint is the solar relative abundance of oxygen to carbon, 
i.e. [C/0] ~0. 



The observed profiles of all doublets can be well de- 
scribed by a model with a single covering factor which takes 
individual values for each ion (smooth curves in Fig. 1). The 
estimated covering factors along with the column densities 
of ions are given in Table 1. 

It is to note that in the present case both the covering 
factors and the column densities are invariant and depend 
neither on the assumed SED nor on the calculation method. 
For instance, a commonly used deconvolution of the line 
profiles into subcomponents delivers the same values for 
the covering factors and column densities as obtained with 
the MCI procedure. 

For the adopted SED (j/~^-^) the ionization curves for 
ions Civ, Ovi and Neviii are plotted by dotted lines in 
Fig. 2 with shadowed area indicating the range of the ion- 
ization parameter U inside the absorber. Both neon and 
oxygen are a-elements and are expected to have the rela- 
tive abundances [Ne/0] ~0. This gives a predicted column 
density A^novih ~ 5 x 10^^ cm~^ with the corresponding 
synthetic profile shown by the dashed curve in the Neviii 
panel in Fig. 1. Obviously the line is too shallow to ac- 
count for the observed Ne viii absorption. This is valid for 
all power law SEDs since the arrangement of C iv, O vi 
and Ne viii ionization curves does not vary with increas- 
ing/decreasing spectral index: the curves are simply shifted 
along the U axis to the right/left, but in the U range en- 
suring the ratio [C/0] ~0 the fraction of Neviii remains 
too small. 

The situation changes for a SED with a break around 
4 Ryd like the one plotted by solid line in Fig. 3. The ion- 
ization curves for this SED are given in the right hand 
part of Fig. 2 (solid lines) with the shadowed area indi- 
cating the range of the fractions of CiV and Ovi ensur- 
ing [C/0] ~0. Now the fraction of Neviii is nearly equal 
to that of Ovi leading to the column density of Neviii 
an order of magnitude higher than from the power law 
SED: iVNeviii ~ 5 X 10^'' cm~^ . The corresponding syn- 
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Table 1. Parameters of the NAL system ^abs = 2.198 towards 
J 2233-606 



Ion 


Ni, cm ^ 


Ci 


Hi 


(1.1 ±0.1)E14 


0.60 ± 0.03 


Civ 


(1.7±0.1)E14 


0.72 ±0.01 


Nv 


(5.3±0.8)E14 


0.70 ±0.02 


Ovi 


(3.4±0.7)E15 


0.82 ± 0.02 




Log U 



Fig. 2. Ionization fractions of different ions calculated for: 
dashed lines - power law ionizing spectrum Fi, oc solid 
lines - ionizing spectrum with a break at ~4 Ryd (Fig. 3). 
Calculations arc performed with solar abundances and for the 
thermal and ionization equilibrium. See text for more details. 



thetic profile for Nevill is plotted in Fig. 1 for the cover- 
ing factor Cncviii = 0.9 (smooth line). A noisy HST/STIS 
spectrum at the positions of Ne VIIIAA770, 780 A prevents 
accurate determination of both the covering factor and col- 
umn density, but it is clearly seen that the adopted - bro- 
ken power law - SED produces A^nc vm consistent with ob- 
servations. This SED gives also the following abundances: 
[C,0,Ne/H] = 0.5, [N/H] = 0.9-1.0, i.e. compared to the so- 
lar values nitrogen is 2.5-3 times overabundant than other 
elements. 

The absence of low-ionization lines and unknown metal- 
licity of the z^bs =2.198 absorber does not allow us to re- 
cover all details of the underlying SED of the ionizing 
radiation. Numerical simulations with different types of 
SED show that the simultaneous existence of significant 
amounts of CiV, Ovi and NeVili can be provided only by 
a sharp break in the spectrum around the energy of 4 Ryd. 
However, other spectral shape features are determined with 
much less accuracy. For instance, positions of the ionization 
curves shown in Fig. 2 are insensitive to the spectral index a 
in the energy range below the discontinuity point A (shown 
in Fig. 3) at 1 Ryd <E <4 Ryd. 

A deeper break after point A does not affect the ioniza- 
tion curve of hydrogen, but shifts the ionization curves of 



ions parallel to the right. The condition [C/0] ~0 is now 
fulfilled at higher U with correspondingly lower fraction of 
neutral hydrogen - thus, we obtain the same column den- 
sity of Neviii but the metal content becomes lower: the 
dashed hue SED in Fig. 3 gives metallicity [C,0,Ne/H] = 
0.18 with stable overabundance of nitrogen by 0.5 dex. 
Power law index between points B and C can vary from 
a = — 0.3toa = 0.1 without significant changes in the pre- 
dicted A^Noviii, but harder spectra deliver more Nevill and 
are probably preferable. The discontinuity point C and the 
power law index beyond it are set in a way to account for 
measurements of AGN/QSO luminosity at E = 2 keV = 
147 Ryd (Steffen et al. 2006): again, metal lines available 
in the Zabs = 2.198 system are weakly sensitive to the coor- 
dinate of this point. 

All calculations above were performed assuming gas in 
thermal equilibrium. However, conditions in the vicinity of 
the central engine vary rapidly, and gas can become over- 
heated. Fig. 4 shows the dependence of the ionization frac- 
tions of different ions on the temperature calculated for 
the UV spectrum from Fig. 3 (solid line) and the ioniza- 
tion parameter U ~ 0.45 which for the solar metallicity 
gives the equilibrium gas temperature of 22000 K. While 
the fractions of CiV and Ovi decrease with rising T, the 
fraction of Ne viii increases, i.e. the overheated gas is more 
abundant in Neviii. On the other hand, the cooling time 
is determined by the gas density, and, hence, the rarefied 
volumes - those traced by Neviii - remain longer in the 
overheated state. Thus, a fraction of Neviii at the value of 
U corresponding to [C/0] ~0 will be higher in the cooling 
gas than that at the equilibrium. 

To summarize, the ionizing spectrum with a sharp break 
in the intensity around 4 Ryd can easily provide conditions 
for significant amounts of C IV, O VI and Ne VIII to arise in 
the same absorbing cloud. 

The results obtained allow us also to make some conclu- 
sions concerning the physical state of the absorbing cloud 
itself. Firstly, the ratio Nciv/^Ovi = 20 along with the 
constraint [C/0] ~0 and the absence (or small amount) 
of O IV (Fig. 1) unambiguously point to photoionization as 
the source of the observed ionization state: in case of coUi- 
sional ionization the leading oxygen ion at the temperature 
which corresponds to Nciv/Nqvi = 20 is Oiv with col- 
umn density an order of magnitude larger than A^ovi (cf., 
Sutherland & Dopita 1993). Previously it was noted by 

The second conclusion concerns the ion-dependent cov- 
ering fac;tors. For the system under study this dependence 
was stated previously by Petitjean & Srianand (1999), and 
they explained this fact by different distances to the broad 
emission line region. However, distribution of covering fac- 
tors as shown in Table. 1 directly follows from a model of 
the associated absorption arising in a gas cloud with fluc- 
tuating density. Namely, ions of higher ionization stages 
trace rarefied gas which can be quite extended whereas 
low ionization ions originate in more dense and, hence, 
compact volumes. This naturally explains the increasing 
value of the covering factor from H I to O VI and the reason 
why the synthetic Ne VIII profile in Fig. 1 was plotted with 

CNeviii = 0.9. 

2.1.2. System at Za.hs = 2.205 

This portion of the associated absorption complex differs 
from the system considered above: profiles of metal doublets 
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Fig. 3. Spectral energy distributions (panel a intensity, panel b 
flux) of the quasar radiation proposed for the associated systems 
towards J 2233-606 to explain the presence of Civ, Ovi, and 
Nevin. The spectra axe normalized so that J^^iKv = 1 Ryd) = 
1. The solid line between 1 Ryd and point A corresponds to 
a = 1.5, the dashed line within the same region - to a = 1.2. 
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Fig. 4. Ionization fractions (in logj^o) of different ions versus ki- 
netic temperature calculated for the incident ionizing spectrum 
shown in Fig. 3 (solid line) at the ionization parameter U = 0.45. 
For the solar metallicity the equilibrium temperature is 22000 K. 



arc more shallow and broad, column densities of C IV and 
O VI seem to be not very different, but the amount of Ne Viii 
is nevertheless significant (Fig. 5). Preliminary calculations 

have shown that profiles of the C ivAA1548, 1550 A linos in 
this system cannot be fitted with a single covering factor. 



Table 2. Column densities and covering factors for the broad 
(b) and four narrow (n) components of the associated system 
^abs = 2.205 towards J 2233-606 shown in Fig. 5 



Ion Ni. VAO.cnr^ C, C, 

Xb) (u) (h) (u) 



Hi <3E13 <0.5 

Civ (2-4)E13 (5.5 - 5.7)E14 0.6 0.2-0.35 

Ovi (1.5-1.7)E15 (1.7-1.9)E15 0.65 0.3-0.4 

Neviii (6-7)E14 0.7 



Winding shape of the profiles indicates that the observed 
absorption may be caused by several gas clouds, each with 
its own covering factor. For a simplest model with multiple 
cloud overlapping the observed intensity can be written as 

n 

Iv = l[(l-Ci,k + Ci,ke-^''^^^) , (2) 
fe=i 

where Ci,fe is a covering factor of cloud k for ion i, and v is 
radial velocity. 

The unsaturated and unblended C IV lines can be de- 
convolved into separate subcomponents. Calculations give 
one broad component with the Doppler parameter 6 = 74 
km s~^ and covering factor Cbroad ~0.6 and four narrow 
components with b =10-26 km and covering factors 
(^narrow =0.2-0.35. Thc column density of the broad com- 
ponent is (2-4) xlO^'' cm^^ , which is only a small fraction 
of a joint column density 5.5 x 10^^ cm~^ of the narrow 
components (Table 2) . Synthetic profiles of the broad and 
narrow components as well as the combined line profile are 
plotted in Fig. 5. Since the profiles of Ovi are saturated 
and smooth they can l)e deconvolved only under additional 
constraints. Assuming the same velocity and component 
structure as for Civ doublet, we obtain the column den- 
sity for the broad Ovi component ^1.5 x 10^''' cm^^ and 
similar column density for the combined narrow compo- 
nents (Table 2). Thus, the broad component shows the ra- 
tio Nqvi/Ncw ^40- 70 which points to a highly ionized 
gas with the ionization parameter U ~1 (Fig. 2), whereas 
thc narrow components have Nqvi/Nqiv ^3 and are much 
lower ionized with U ~0.1. With these parameters the pre- 
dicted column density for NeVlll is ~6 x 10^^ cm~^ for the 
broad component (plotted in Fig. 5 for the covering factor 
0.7) and ^10^'^ cni^'^ for the narrow components. 

The hue N vA1238 A is blended with CaiiA3968 A from 
^abs ~0 which hampers its deconvolution into components 
and estimation of their covering factors. However, trial cal- 
culations with different constraints show that an overabun- 
dance of nitrogen is present in both broad and narrow com- 
ponents. 

2.2. NAL systems towards HE 1341-1020 

The spectrum of the quasar HE 1341-1020 was obtained 
with the UVES/VLT in the framework of the ESQ Large 
Program 'QSO Absorption Line Systems' (ID No. 166. A- 
0106). Data reduction was performed by B. Aracil. 

HE 1341 -1020 is a mini-BAL quasar with broad absorp- 
tions in Civ, Nvi and Ovi extending over 3000 km s~^ . 
There is also a narrow absorption system at 2;abs = 2.147. 
The quasar emission redshift estimated from the C iiA1334 
A and C IIIA1990] A emission lines is ^em = 2.1485 ±0.002. 
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Fig. 5. Same as Fig. 1 but for the Zabs = 2.205 system towards J 2233-606 (sohd-hne histograms). The zero radial velocity is fixed at 
z = 2.2050. Synthetic profiles of C iv and O vi are produced by broad and shallow absorption (dashed curves) with covering factor 
C ~0.6, and by several overlapped narrow absorptions (dotted curves) with C ~0.3. Smooth curves represent the convolution of all 
components according to Eq.(2). The dashed line in the Neviii panel - predicted profile for the broad component corresponding 
to the ionizing spectrum with a break at ~4 Ryd (Fig. 3, solid line). The dashed line in the Ly-a panel - predicted upper limit 
for the broad component. 



2.2.1. System at Zabs = 2.147 

Ionization state and restored SED. The system at Zabs = 
2.147 exhibits a wealth of lines of metal ions ranging from 
low ionization Si , N II, Mg II up to high ionization N V 
and Ovi (Fig. 6). The presence of many elements in dif- 
ferent stages of ionization makes it possible to reconstruct 
the shape of the HE 1341-1020 spectrum in more detail 
than that previously described for J 2233-606. Trial calcu- 
lations revealed that power law spectra as well as the AGN 
spectrum of Mathews & Ferland (1987) failed to reproduce 
the observed ratio Nii/Nv resulting in the synthetic line 
profiles which significantly underestimate the observed in- 
tensity of N IIA1083 A, and yielded relative element abun- 
dances which did not comply with existing observational 
and theoretical data (e.g., [Mg/0] ~ 0.7, [Si/0] ~ 0.3). 
Several types of SEDs with the break at 4 Ryd were tried 
as well and as a result the SED recovered from the Zabs — 
1.78 associated system towards HE 0141-3432 (Reimers et 
al. 2005, Sect. 4.2) was chosen as the best initial guess for 
the spectral shape recovering procedure. The search for an 
optimal spectral shape was directed by the requirement to 
reproduce all observed lines under the constraint [Si ,Mg , 
0/C] < 0.2. Details of calculations are given in Reimers 
et al. (2005, appendix). 

The resulting SED is shown in Fig. 7, and the corre- 
sponding synthetic profiles are plotted by solid curves in 
Fig. 6. The estimated physical parameters and column den- 
sities are given in Table 3. 

As it was for the previous system, the spectral shape 
again has a sharp break at ~4 Ryd, but in the present case 
it is also possible to fix the depth of this break and the 
broken power law slope between 1 and 4 Ryd. 

The column density of neutral hydrogen is A^hi ^ 10"'^'' 
cm~^ , i.e. the ^abs = 2.147 system is a sub-LLS (optical 
depth in the Hi Lyman contimmm r^' < 0.6). The pre- 
dicted column density for Hen is Ajjen ^ 3 x 10^® cm~^ 
which makes the system moderately optically thick to the 
Hell Lyman continuum (r^''" — 4.7). 

The ionizing spectrum is reconstructed from the ob- 
served absorption lines and therefore is sensitive to the local 
processes in the absorbing cloud itself (for more details, see 



Sect. 3 in Agafonova et al. 2005). In particular, for the case 
in question the restored quasar spectrum can be affected by 
the He li Lyman continuum absorption inside the cloud, i.e. 
the incident SED can be harder (probably by 0.2-0.3 dex) 
at energies E > A Ryd than the spectrum shown in Fig. 7. 
The He ii Lyman continuum absorption is accompanied by 
the recombination emission line of He ii Ly-a (304 A) and 
by the two-photon emission which both produce a charac- 
teristic emission feature at S < 3 Ryd in the transmitted 
spectrum of the incident quasar radiation. Whether such a 
feature is present or not in the restored spectrum cannot 
be stated unambiguously because of the saturated profiles 
of the Si III and Si IV lines (fractions of silicon ions are most 
sensitive to the spectral shape of the ionizing radiation at 
E <3 Ryd as seen in Fig. 2 in Agafonova et al. 2007) . 

The measured mctallicity of the absorbing gas is solar 
(carbon) or slightly oversolar (a-elements) with clear un- 
derabundances of nitrogen (0.6 solar) and aluminium (0.3 
solar). 

Distance to the QSO. The presence of the ground fine- 
structure line Cii*A1335 A arising from the first excited 
•^P^ level provides a rare opportunity to measure from the 
photoionization model the mean gas number density and, 
hence, the distance to the light source. For an ion in the 
interstellar (intcrgalactic) medium, the ratio of excited (^2) 
to ground-state (ni) population is equal to the ratio of the 
collisional excitation rate (3i_>2 to the spontaneous transi- 
tion probability ^2^1 (Bahcall & Wolf 1968): 

- = (3) 

The atomic data for C 11* are the following (Silva & Viegas 
2002): Aa^i = 2.291 x lO"*' s-\ the excitation rate by 
coUisions with electrons at Tkin = 10** K ql^2 — 1 x 10"'^ 
cm-^ s~^. Since collisions with other particles have much 
lower excitation rates, we put Qi^2 = 9i^2 '^e- 

Prom the measured column densities A^cn and A^cn* 
(see Table 3) we obtain the ratio of the mean values n2 to 
m 
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Fig. 6. Same as Fig. 1 but for the Zabs = 2.147 system towards HE 1341-1020 (sohd-line histograms). The zero radial velocity is 
fixed at z = 2.14736. Synthetic profiles corresponding to the ionizing spectrum with a break at ~4 Ryd (Fig. 7) are plotted by 
smooth curves. Note different vertical scales in some panels. 



Since this system has a high degree of ionization 
{Nui/Nh = n^+/nii 1, see Table 3), and collisions with 
other particles have much lower excitation rates, the total 
gas density equals to no — 0.5 cm~^ (the contribution of 
the ionized helium is ignored since it has a small effect). 

If this system would be an intcrgalactic absorber, then 
the ionization parameter U = 0.02 would correspond to the 
gas number density no ~ 0.001 cm~^ (assuming the inten- 
sity J912 of the metagalactic UV background radiation at 
z = 2.0 as given by Haardt fc Madau, 1996). This is a clear 
indication that the ionizing background at the position of 
the Zabs = 2.147 system is enhanced by more than 2 orders 
of magnitude as compared to the intergalactic background. 

In order to estimate the distance to the light source the 
QSO continuum luminosity at the hydrogen Lyman limit 
is to be known. The intrinsic luminosity £,y^ can be 
determined from the comparison of 1^ = 17.1 with the spe- 
cific flux of a star having my = 0.0 outside the Earth's 



atmosphere, Fy ~ 3.8 x 10 ^" ergs cm ^ s ^ Hz ^ : 

/v = 10"° '' "-^F^ = 5 X 10^2^ ergs cm'^ s'^ Hz'^ 

With the approximation for luminosity distance c^l pro- 
posed in Riess et al. (2004), one obtains c^l — 13 Gpc and 
the apparent monochromatic luminosity 

£^ = 47r 4 /v (1 + z)-^ ~ 3 X 10^^ ergs s"^ Hz"^ , 

which results in the luminosity near the Lyman limit 

= ' :^ 1.5 X 10^1 ergs s'^ Hz'^ . 

Given the gas number density no, the distance r from the 
QSO to the absorbing cloud can be calculated from the 
ionization parameter J/q which is defined as 

„. = J3S!L = >i, ,4, 

47r c no no 
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Fig. 7. SED estimated from the 2:abs = 2.147 system towards the 
quasar HE 1341-1020 (zem = 2.1485). 



according to Eq.(2) gives for the broad component the col- 
umn density iVciv,broad = 4.2 x 10^^ cm~^ and the covering 
factor Cc IV, broad = 0.8, and for the four narrow compo- 
nents the total column density A'^c iv, narrow = 5-8 x 10"'^^ 
cm~^ and the covering factor 0.16 < Cciv, narrow ^ 0.36 
(thin smooth and dotted curves in panel CiV in Fig. 8). 

From the doublets N v and O vi only the hues N vA1238 
A and OviAlOSl A are available. Assuming for Nv sim- 
ilar velocity structure and covering factors as calculated 
for Civ, one obtains for the broad Nv component the 
column density A^nv, broad = (1-4 — 1.6) x 10^^ cm~^ . 
The saturated OviAlOSl A line indicates the covering fac- 
tor of Covi.broad — 0.90 — 0.95, which gives the upper 
limit for the column density of the broad O VI component 

A^Ovi.broad < 2 X lO^^ ^m-Z 

The ratio A^o vi,broad/-/Vciv,broad ^ 50 together with the 
assumption of solar relative abundances [C/0] ^0 lead to 
the range of the ionization parameters U which translate 
the column densities of C iv and N v observed in the broad 
component into the reliable relative overabundance of nitro- 
gen [N/C] > for both a power law SED and a SED with a 
discontinuity at 4 Ryd. Thus, the gas comprising the broad 
component in the mini-BAL absorption at the radial veloc- 
ity V ~ —4000 km s~^ obviously differs from the gas pro- 
ducing the narrow line absorption at Zabs = 2.14736 (Fig. 6) 
where the nitrogen under abundance of [N/C] = —0.15 was 
detected. This will be discussed further in Sect. 3.2. 



where 
Q(H°) 



- — av 
hv 



(5) 



is the number of hydrogen ionizing photons emitted per 
unit time by the central source. 

With the estimated above Lyman continuum luminos- 
ity one finds Q(^^) ~ 2 x 10"''^ photons s^^, assuming 
Ci, = {v/vc)~^- A substitution of the numerical values 
in (4) provides r ~ 240 kpc. Accounting for the optical 
depth r < 0.6 at the hydrogen ionization edge due to ab- 
sorption within the cloud, the distance can be reduced to 
~170 kpc. 

It is unclear whether an absorber at such a large dis- 
tance from the light source can be gas ejected from the 
QSO host galaxy. However, distances as large as hundreds 
of kiloparsecs (mostly estimated as lower limits) arc not 
exceptional for the associated systems (Morris et al. 1986; 
Tripp et al. 1996; D'Odorico et al. 2004; Reimers et al. 
2005). 



2.2.2. System at Zabs = 2.107 

Mini-BAL absorptions in Civ, Nv and Ovi at -7000 
km s~^ < V < —1000 km s~^ are formed by a large number 

of overlapping components (Fig. 8) which in most cases do 
not allow to estimate covering factors and, hence, to mea- 
sure accurately column densities. Fortunately, portions of 
the CivA1548 A profile between -4100 km s^^ and -3855 
km s^^ (zabs — 2.1066) and, correspondingly, CivA1550 
A profile at —3638 km s~^ < u < —3396 km s~^are not 
blended and can be used for quantitative analysis. A wind- 
ing shape of the line profiles favors the model used above for 
the Zabs = 2.205 system towards J 2233-606; several narrow 
components seen against the broad one. Component fitting 



2.3. NAL systems towards HE 2347-4243 

An associated complex towards the quasar HE 2347-4243 
(zom = 2.902) consisting of many absorption systems and 
extending over 1500 km s^^was described in detail in 
Fechner et al. (2004). It was found that column densities 
of ions measured in certain systems favored the ionizing 
spectrum with a break at E ^4 Ryd. However, an exact 
spectral shape in the range E > I Ryd was not estimated. 
Below we re-analyze a few systems from this complex with 
the objective to reach more tenable conclusions about the 
SED of HE 2347-4243. 



2.3.1. System at z^hs^ 2.898 

This system (component 8 in Fechner et al. 2004) exhibits 
strong lines of the doublets Civ, Nv and Ovi which allow 
us to estimate accurately the covering factors and column 
densities (Fig. 9). There is a weak hydrogen absorption Ly- 
a, whereas the observed intensity at the position of Ly-/3 
may be due to a blend with some forest line. However, these 
lines can be used to set a limit on the covering factor for 
Hi. The CiiiA977 A fine is also observed, albeit blended. 

The measured column densities and covering factors for 
ions Civ, Nv and Ovi are given in Table 4. The appar- 
ent intensity at the position of Ly-/3 constrains the cov- 
ering factor and column density for Hi: Cm < 0.6, and 
A^Hi < 2.3 x 10^^ cm~^ . The covering factor for Cm can 
lie between Cm and Cciv, i-e. 0.60 < Ccm < 0.97. The 
observed intensity at the expected position of the C ill 
line allows us to estimate the upper limits to its column 
density: A^cm < 9.0 x 10^^ cm'^ if Ccm = 0.60, and 
Ncui < 4.7 X 10^2 cm-2if Cc„, = 0.97. 

Taking the ratio Nqiv/Nqvi = 2.3 and assuming so- 
lar relative abundance of carbon to oxygen, [C/0] ~0, we 
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Fig. 8. Mini-BAL absorptions towards HE 1341-1020. Two upper panels represent portions of the original QSO spectrum in the 
vicinity of broad emission lines H I Ly-a and O Vi+H I Ly-/3. The smooth curves here are the local continua calculated through the 
spline interpolation between 'clear' continuum windows seen at high spectral resolution (FWHM ~6 km s~^ ) in the UVES/VLT 
data. The lower panels show the zoomed portions marked by long brackets in the upper panels. The vertical axis is normalized 
intensities. The zero radial velocity is fixed at Zcm— 2.1485. Two combined panels with Hi A1216 profiles are, respectively, the 
general view with labeled absorption lines identified at different redshifts and the same but zoomed in the vertical scale spectrum 
to illustrate broad and shallow absorption features (marked by numbers over horizontal lines) caused by Hi Ly-a arising from 
a fast gas outflow. The corresponding broad and shallow features are seen in the O vi panel. Absorption centered aX v = —4000 
km s~^(zabB= 2.107) can be deconvolved into one broad and shallow component and several overlapped narrow components 
(shown by dotted curves). Smooth curves show: CiV - synthetic profile of the convolved broad and narrow absorptions; Nv - 
broad component corresponding to the covering factors and the velocity structure as determined from C iv; O VI - limiting profile 
for the broad component. Long brackets in the three bottom panels mark components of the corresponding doublets to illustrate 
their strong overlapping. 



can estimate the ionization parameter U for a given SED. 
The limits set above on A^cm reject unambiguously any 
type of a power law spectrum as well as the AGN spec- 
trum of Mathews & Ferland (1987): for U corresponding to 
[C/0] --^O these spectra significantly overpredict Ncm- 

In order to comply with the limits on Ncm the ion- 
izing spectrum should be quite hard at 1 < < 4 Ryd 



{a < 1.0) and have a sharp break at E ~4 Ryd with the 
intensity depression by more than an order of magnitude. 
The spectrum of HE 2347-4243 is known up to the wave- 
length of about 350 A (2.6 Ryd) and it indeed shows a very 
hard continuum in the EUV range with a = 0.56 (Rcimcrs 
et al. 1998; Fechner et al. 2004). Adopting this index, we 
can estimate the depth of the intensity break at 4 Ryd: 
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Ja/Jb = 20 (assuming Ccm = 0.60), and Ja/Jb = 30 
(assuming Ccm = 0.97). 

With these EUV spectra, the carbon and oxygen abun- 
dances in the Zabs = 2.898 system can vary between 5-15 
solar - depending on the adopted covering lactor for H I 
and the depth of the continuum break at 4 Ryd, but nitro- 
gen remains always under abundant: [N/0,C] < 0. Absolute 
values of the abundances are uncertain not only due to un- 
known covering factor of Hi, but also because of a possible 
overionization of Hi as a result of high temperature: calcu- 
lations with CLOUDY shows that the H i fraction decreases 
more rapidly than the fractions of CiV and Ovi when 
temperature rises above the thermal equilibrium level. On 
the other hand, time to reach the ionization equilibrium 
is longer for Hi than for other ions (Osterbrock 1974). It 
means that calculations which use the Hi fraction corre- 
sponding to the thermal cquilibriTim can significantly un- 
derpredict the total hydrogen amount and deliver artifi- 
cially boosted metallicity. In principle, a neighbor system 
centered at —415 km s~^ (Fig- 9) may be considered as a 
support for this statement: although accurate calculations 
are prevented by line blending, an upper limit for the gas 
metallicity does not exceed 3 solar values. Unfortunately, it 
is impossible to distinguish between two interpretations of 
absorption lines seen at t; < —325 km s~^ : whether they are 
caused by the overlapping of clouds with different metallic- 
ities or by clouds with similar metallicity but some of them 
not in thermal equilibrium. 

2.3.2. System at Za.bs = 2.901 

This system (centered at w ~ —85 km s~^ in Fig. 9, com- 
ponent 10 in Fechner et al. 2004), exhibits unblended weak 
lines of the doublets Civ and Nv along with strong sat- 
urated lines of Ovi partially blended in the wings. The 
line Hi Ly-a is saturated and black at the center which 
means that the absorber completely covers the light source. 
Assuming constant metallicity across the cloud, it is pos- 
sible to reconstruct the line profiles from their unblended 
parts. The obtained column densities are given in Table 4. 
A spectral feature at v ~ —25 km s~^ in the Ovi profile 
is due to the assumption of the constant metallicity and is 
uncertain since the corresponding parts in both O VI lines 
are blended. However, its input to the total column density 
of O VI is insignificant because A'^o vi is mostly determined 
by the central portion of the Ovi lines which are clear. 

The measured ratio Nqvi/^Civ = 200 is extremely high 
and points to the ionizing radiation much stronger than the 
intergalactic: for the intergalactic absorbers with similar 
column densities of Hi the ratio A^ovi/-^Civ is about 10- 
13 (Reimers et al. 2006). Thus, the system under study is 
probably located close to the quasar. 

With the same ionizing spectrum as was determined for 
the Zabs = 2.898 system, the ionization parameter U can 
be estimated if the relative abundance of carbon to oxygen 
[C/0] is known. It is interesting to note that the condi- 
tion [C/0] ~0 used to fix U in the preceding sub-sections, 
cannot be fulfilled in the present case: starting at U ^1 the 
ionization curves of O VI and C iv are almost parallel giving 
the ratio of the ion fractions lg(Tovi/T"civ) = 0.8 and thus 
producing a stable overabundance of oxygen [0/C] ~0.2. At 
lower U { <1) the overabundance of oxygen increases. 

From measurements in the Galactic and cxtragalactic 
H II regions it is known that low-metallicity gas can be ov- 



Table 3. Physical parameters of the Zabs ~ 2.147 metal absorber 
towards HE 1341-1020 derived by the MCI procedure (uncer- 
tainties of the physical parameters are 10-15%: Uo - mean ionizar 
tion parameter; Nh - total hydrogen column density; crv,Cy - 
dispersions of velocity and density distributions, respectively; 
Zx = ATx /Nu - abundances of the individual elements) 



Parameter 



"[Zx] = log(Arx/iVH) - log(iVx/iVH)0 . 
'Calculated using the velocity and density distributions 
derived from hydrogen and metal profiles. 



erenriched in oxygen with the safe upper bound [0/C] < 0.5 
(Henry et al. 2000; Nava et al. 2007). This gives U > 0.2 
and metal content relative to solar [C] < —1.0, [O] < —0.5, 
[N] < —1.2. These relative abundances may indicate SNe II 
explosions as a main source of metal enrichment in the 
Zabs = 2.901 absorber. 

Metal content of the absorbing gas redward to the Zabs = 
2.901 absorber (systems centered at v = km s~^ and 
V = 115 km s~^ in Fig. 9) is uncertain due to blending 
of both Ovi lines with forest absorptions. However, the 
column density of neutral hydrogen can be estimated quite 
accurately from the break at 912 A (rest frame) clearly 
seen in the spectrum of HE 2347-4243: N{Rl) = (1.6 - 
2) X 10^^ cm~^ . Depending on the ionization parameter 
of the gas and the spectral shape of the incident ionizing 
radiation the corresponding amount of He ii ranges in the 
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Table 4. Column densities and covering factors for the associ- 
ated systems at Zabs = 2.898 and 2.901 towards HE 2347-4243 
shown in Fig. 9 





z = 2.898 


z = 2.901 




Ion 


N, cm-2 C 


N, cm-2 


C 



Hi 


2.3E13 
1.1E13 


0.6 
0.97 


(7.6 ± 1.0)E14 


1.0 


Cm 


9.0E12 
4.8E12 


0.6 
0.97 






Civ 


(1.2±0.1)E14 


0.97 


(6.7±1.0)E12 


1.0 


Nv 


(2.6±0.2)E13 


0.97 


(1.5±0.3)E13 
(1.3±0.2)E15" 


1.0 


Ovi 


(2.8 ± 0.2)E14 


0.985 


1.0 



"Both O VIAA1031, 1037 A lines are blended in the red wings 
which are restored from v-d distributions obtained from 
other metal lines and assuming constant metallicity within 
the absorber. The column density is set by the portion 
—60 < V < 40 km s~^ which is unblended. 



Table 5. Column densities and covering factors for the associ- 
ated system at Zabs = 2.352 towards Q 0329-385 shown in Fig. 11 



lou 



A'. 



Hi 


(1.3 - 1.4)E13" 


0.6 




(7 - 8)E12'' 
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Cm 


3.6E12'' 


1.0 


Civ 


(3.2±0.1)E13 

2.9E12'' 


1.0 


Niii 


1.0 


Nv 


(4.0 ± 0.2)E13 


1.0 


Ovi 


(1.2 ± 0.1)E14 


1.0 



"The red wing of the apparent profile inconsistent with 
the assumption of constant metallicity across the 
absorber. The synthetic profile is restored on base of 
velocity-density (v-d) distributions obtained from metal 
lines and assuming homogeneous metallicity. 
''Estimated from v-d distributions obtained from the 
metal lines C IV, N v, O Vi and SED shown in Fig. 10. 



interval (2 — 5) x lO^'^ cm~^ and, thus, Hell can soften the 
transmitted UV spectrum by 0.3-0.4 dcx at i? > 4 Ryd ^. 
Since relative location of the absorbers along the line of 
sight from the quasar towards the observer is unknown, it 
cannot be excluded that the absorber at 2;abs = 2.898 follows 
in physical space after the absorbers considered here. This 
means that the SED reconstructed from the Zabs = 2.898 
system may represent the quasar radiation softened partly 
by the intervening gas seen at —150 km s~^ < t; < 150 
km s~^ . 

2.4. System at ^abs = 2.352 towards Q0329-385 

The spectrum of the quasar Q 0329-385 was obtained 
with the UVES/VLT in the framework of the ESO Large 
Program 'QSO Absorption Line Systems' (ID No. 166. A- 
0106). Data reduction was performed by B. Aracil. The 
system at Zabs = 2.352 was firstly described in Bergeron et 
al. (2002) in their study of intergalactic O vi absorbers and 
recently by Schaye et al. (2007). 

The redshift of the quasar Q 0329-385 is z^m = 2.435 
(Hi Ly-a) which means that the Za,bs= 2.352 system is 
detached from the QSO by ^^7400 km s"'^ . There are clear 
lines of the doublet Civ, whereas doublets Nv and Ovi 
are blended each in one component (Fig. 11). The relative 
intensities of both C IV lines do not show any evidence of the 
incomplete coverage of the light source. The system is very 
much ahke to that at Zabs = 2.898 towards HE 2347 4243 
(Sect. 2.3.1): strong lines of highly ionized metals combined 
with a weak Hi line. 

The measured column densities are given in Table 5. 
The present system shows the column density ratio 
Nqvi/Nqiv = 3.7 which is only slightly higher than 
■^Ovi/-^Civ = 2.3 of the 2abs= 2.898 system, but it con- 
tains much more nitrogen: N^y/Nciv = 1-34 compared to 
-/VNv/iVciv = 0.21 at Zabs= 2.898. 

The apparent profile of Ly-a is inconsistent with the as- 
sumption of constant metallicity throughout the absorber 
(synthetic profiles are shown by the smooth curves in 

1 Smette et al. (2002) suggest for this system Ar(Hell) > 2.4 x 
10^* cm^^ , which is a clear ovcrcstimation: the extremely hard 
at 1 < -E < 4 Ryd spectrum of He 2347-4243 cannot provide 
ri = iV(Heii)/iV(Hi) > 100 for any value of the break at 4 Ryd. 



Fig. 11)^. The Hi Ly-/3 line is blended and can be em- 
ployed merely to set a lower limit for the covering factor 
and, correspondingly, an upper limit for iV(Hl). Since the 
ClllA977 A line is blended with a strong Hne CivA1548 
A from the Za,hs= 1-115 system, the only way to fix the 
ionization parameter is to use the condition [C/0] ~0. 

With this restriction, any reasonable spectral index a 
of a pure power law SED {a = 1.0 — 1.8) leads to the 
C and O abundances of 10-20 solar and to the nitrogen 
abundance of 20-40 solar. These values seem to be definitely 
too high, even if one assumes that hydrogen is overheated 
(in fact, there are no signs of possible overheating in this 
system). On the other hand, for the ionizing spectra with 
a break at 4 Ryd the condition [C/0] ~0 is fulfilled at 
higher U than for the power law spectra, resulting in a lower 
fraction of neutral hydrogen and, hence, a higher amount 
of total hydrogen whereas the fractions of C IV and O vi 
do not change. Thus, the metal abundances decrease: the 
softer the spectrum at E > 4 Ryd the lower metallicity. For 
example, the spectrum restored for the ^abs = 2.898 system 
(Fig. 10, soHd Hne) gives [C] = [O] = 0.5 and [N] = 0.9. 
The accurate value of the continuum depression at 4 Ryd 
cannot be estimated since the real metallicity of the gas 
cloud at ^abs = 2.352 is of course unknown. However, we 
can conclude that the SED with a break at 4 Ryd is clearly 
preferable for this system. Another conclusion is a revealed 
significant relative overabundance of nitrogen, [N/C] ~0.3- 
0.4, which was constantly reproduced for any tried ionizing 
spectrum. 

3. Discussion 

3.1. The origin of the break at 4 Ryd 

All NAL systems considered above as well as NAL systems 
towards HE 0141 3932 (Reimers ct al. 2005) are best de- 
scribed with ionizing spectra having a break at £^ « 4 Ryd. 
It seems natural to relate this break to the Lyman contin- 
uum absorption in Hell. Then the question arises where 
does this absorption occur? 

^ The velocity shift between the dominant component in H I 
and metal lines was noticed for the Zaha = 2.352 system also by 
Shaye et al. (2007). 
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Fig. 9. Same as Fig. 1 but for the «abs = 2.898 {v ~ —317 km s~^ ) and Zabs = 2.901 {v ~ —85 km s""' ) systems towards HE 
2347-4243 (solid-line histograms). The zero radial velocity is fixed at z = 2.9021. Synthetic profiles corresponding to the ionizing 
spectrum with a break at ~4 Ryd (Fig. 10) are plotted by smooth curves. See text for details. 



Firstly, this absorption can be produced by the NAL 
systems themselves by those with Af(Hi) > 10^® cm~^ . 
This possibility was suggested in Reimers et al. (1997) and 
then discussed also in Smette et al. (2002) and ShuU et 
al. (2004). Among considered absorbers, to this group be- 
long the system at Zabs = 1.7103 towards HE 0141 3932 
(Reimers et al. 2005), and systems at 2;abs = 2.901 (HE 
2347-4243) and ^abs= 2.147 (HE 0141-3932) from the 
present work. The Zabs = 2.147 system with iV(Hl) « lO^'^ 
cm~^ and A''(He ii) w 3 x 10^* cm~^ can noticeable soften 
the incident radiation due to internal absorption in the 
He II Lyman continuum, but the metal lines detected show 
that the incident radiation should already have a signifi- 
cant break a,t E = 4 Ryd (Sect. 2.2.1). The Zabs = 1.7103 
and Zabs= 2.901 systems have 7V(Heii) < 5 x 10^'' cm~^ 
and also cannot account for the whole depth of the break 
at 4 Ryd in the reconstructed spectra. Besides, lines of 
sight towards J 2233-606 and Q 0329-385 with the as- 
sociated systems at Zabs = 2.198 (Sect. 2.1.1) and Zabs = 
2.352 (Sect. 2.4), both strongly favoring the ionizing spec- 
trum broken at 4 Ryd, do not show any absorption with 
7V(Hi) - 10^^ cm-2 in the range z^m - z < 0.3. Thus, 
in general the i;abs ~ -^em systems can hardly be the main 
source of the observed He II Lyman continuum absorption. 

More probably, the bulk of the 4 Ryd discontinuity is 
caused by the He ll absorption in quasar winds which show 
a common occurrence both in quasars and AGNs. The ac- 
cretion disk winds arc supposed to be accelerated by the 
combination of radiation-pressure and magneto-centrifugal 



forces and comprise dense gas near the equator seen as UV 
broad absorption lines (BAL) and highly ionized shielding 
gas which protects the inner wind regions from overion- 
ization by X-rays from the central source (Murray et al. 
1995; Proga et al. 2000; Everett 2005). All quasars consid- 
ered here do not belong to BAL quasars (only HE 0141 
3932 can be classified as mini-BAL). The natural question 
is whether their spectra contain any features which allow 
to conclude that the wind absorption is indeed present and 
that the column density of He ll in the outflowing gas is of 
order -^10^^ cm~^ (the depth of the 4 Ryd break ~1 dex). 

Model calculations with CLOUDY show that for any 
incident ionizing spectrum the ratio r] = A^(Hell)/iV(Hl) 
firstly increases with the rise of the ionization parameter 
U, but attains the constant value when U becomes large 
(Fig. 19 in Agafonova et al. 2007). For power law spec- 
tra this ratio equals 20 — 40 depending on the slope of the 
spectrum (the steeper spectrum the larger 77) and remains 
almost constant at J7 > 0.1. Thus, TV (He 11) of ~10^^ cm~2 
corresponds to (H i) of a few times 10^^ cm~^ . The ion- 
ization parameter U can be constrained if we assume the 
total column density of hydrogen A''(H ) ~10^^ — 10^^ cm~^ 
which is estimated from X-ray observations of quasars 
(Murray et al. 1995; PiconceUi et al. 2005). This gives U ~1- 
5. In this U range, among ions observable in the UV the 
highest fractions exhibit Ovi (1/30 - 1/100 of all oxygen) 
and Nevi, Nevii, Neviii (1/4 -1/10 of all neon) whereas 
the ionization fraction of C iv is negligible. Gas with such 
parameters is expected to leave in a UV quasar spectrum 
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Fig. 10. SED restored from the Zabs ~ 2.898 system towards 
the quasar HE 2347-4243 (zem = 2.902). The dashed hne repre- 
sents the ionizing spectrum corresponding to the covering factor 
Cciii = 0.6, solid Hne - to Com = 0.97 (see text for details). 



its footprints seen as Hi Ly-a A1215 A absorption, Hi 
break at A < 912 A and - in case of high metaUicity - 
Ovi AA1032,1037 A, Nevi AA401,430 A Nevii A465 A, 
and Ne viii AA770, 780 A absorption features. Since the gas 
velocity range in the accretion disk winds exceeds 10000 
km s~^ , these absorption features should look like wide 
and shallow troughs, and a break at the hydrogen ioniza- 
tion edge - as a slow continuous depression towards higher 
frequencies. It is to note that in the given U range the oxy- 
gen exists mostly in the form of Ovii (almost 90% of all 
oxygen) and O viii (~10% of oxygen), i.e. is best accessible 
in the X-ray region. 

We do find the expected features in the quasar spectra 
considered in the present work. The depression at A < 912 
A (rest-frame) is clearly seen in the flux-calibrated spec- 
trum of J 2233-606 shown in Fig. 1 in Sealy et al. (1998) 
and its depth (r^' ~0.3) corresponds to TV (Hi) ~4 x 10^^ 
cm~^ . Depression occurring in the flux-calibrated spectrum 
of HE 2347-4243 (Fig. 3 in Fechner et al. 2004) accounts 
for A/'(Hl) from 3 to 7 x 10^^' cm^^ . For HE 1341 1020, the 
corresponding wavelength range is beyond the observable 
frame, whereas for Q 0329-325 the flux-calibrated spectrum 
is not available. It should be noted that the flux-calibration 
is the mandatory requirement to detect this slow rolling 
continuum depression because otherwise it is disguised by 
instrumental effects. Further on, shallow troughs of hydro- 
gen Ly-a and Ovi absorptions extending over thousands 
km are present in spectra of HE 0141-1020 (Fig. 8, 
marked by horizontal hues), HE 2347-4243 (Fig. 12) and 
Q 0329-385 (Fig. 13). In spite of very low optical depth of 
Hi Ly-a (r^^ ~0.03), the corresponding portions can be 
clearly distinguished in the spectra due to high signal-to- 
noise ratio and high spectral resolution which allows us to 
trace 'clear continuum windows' between absorption fea- 



tures over a wide spectral range. We also checked whether 
the detected troughs could be due to mismatch of echelle 
orders. The answer is no: in the spectral range considered 
the orders are small (~300 km ), and a depression of the 
width of ~1000s km s^^ covers several subsequent orders. 

The optical depth of an absorption line arising in a 
gas flow moving with a large velocity gradient is given by 
(Sobolev 1947, 1960; Castor et al. 1975): 



r = 2.654 X 10'^^ fXn/ 



dv 
ds 



(6) 



where / is the oscillator strength for absorption, A the cen- 
tral wavelength (in A), n the number density (in cm""^ ) 
of atoms/ions in the ground state, v the wind speed (in 
km s~^ ), and s the distance (in cm) along the line of sight. 
The population of the upper excited level is ignored. 

For Hi Ly-a and Ovi A1031 (as an example of a heavy 
element) the corresponding optical depths are 



= 1.34 X 10-12 Thi/ 



dv 
ds 



and 



T^osl = 3.63 X 10-13 riH Tovi/ 



dv 
ds 



where Zq is the oxygen abundance and Thi, Tovi are the 
ionization fractions of Hi and Ovi, respectively. 

The viable estimation for the velocity acceleration dv/ds 
at V > 1000 km s^^ is 10"^ - 10"® km s-^ cm-^ (e.g., 
Murray et al. 1995) with lower acceleration corresponding 
to higher velocities. Thus, to produce the optical depth 
of ~0.03 the number density of neutral hydrogen should 
be riHi ~(2-20)xl02 cm-"^ . For an ionization parameter 
U r^l-b this gives the gas number density of nn = 10* — 10^ 
cm-3 . Absorption in a gas with lower density will be indis- 
tinguishable from noise fluctuations in the quasar contin- 
uum windows whereas more dense parts of the disk wind 
will produce pronounced absorption which can be seen not 
only in H I, O VI, and Ne viii but - due to lower i7 - in such 
ions as C iv, N v or even Si iv as well. Large density fluctua- 
tions are expected in the disk wind due to flow stratification 
and due to the presence of various instabilities which give 
rise to shocks and dense shells as they propagate from the 
wind base outward and grow rapidly to non-linear regime 
(Owocki et al. 1988: Proga et al. 2000). Thus, in general, 
the wind absorption is represented by the structure stretch- 
ing over 1000s km s^^ and consisting of the low-contrast 
continuum depression intermitted by parts without any ab- 
sorption and by deeper absorption features caused by both 
density and velocity fluctuations'^. In the quasar spectra 
considered here the broad and shallow absorption features 
can be detected at radial velocities up to 20000 km 
from the quasar's systemic redshifts which corresponds to 
Az ~0.2 at 2 < 2; < 3. Thus, the superposition of ab- 
sorption lines from the intervening Ly-a forest clouds is 
inevitable. Absorptions of Hi and Ovi in Fig. 8, 12, 13 
show all elements of this picture*. 



^ cf. Levshakov & Kegel (1998), where a 'line-like' struc- 
ture caused by velocity fluctuations superposed on the general 
Hubble flow was calculated in the approximation of the constant 
gas density and neglecting the coupling of the velocity field with 
the radiation field. 

Similar complex absorption-line system consisting of shallow 
and broad {FWHM -700 km s"^ ) and narrow {FWHM < 20 
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Fig. 11. Same as Fig. 1 but for the Zabs = 2.352 system towards Q 0329-385 (sohd-line histograms). The zero radial velocity is 
fixed at z = 2.35205. Synthetic profiles are plotted by smooth curves. 



Another important question is to what extent the 
quasar's outcoming radiation is affected by the Hell Lyman 
continuum absorption. It seems quite probable that outflow 
zones which produce BALs arc opaque in He ii continuum. 
X-ray observations of BAL quasars show that X-ray ab- 
sorbing gas has A''(H ) = 10^^ - 10^^ cm'^ (Gallagher et 
al. 2002). The BAL systems are mostly identified by strong 
C IV absorption which supposes that the ionization param- 
eter U is less than 1 for a power law incident continuum 
(or that incident continuum is soft at E 4 Ryd). This 
gives the column density of neutral hydrogen A''(H)~10^'' 
cm~^ and of single ionized helium A'^(He ii) ^ few times of 
10^* cm^^ . Thus, the quasar spectrum transmitted through 
the BAL-producing region should have a pronounced break 
(more than 1 dex) at 4 Ryd. However, BAL quasars com- 
prise only about 15% of quasars population (Reichard et 
al. 2003) what means if one considers the BAL effect as 
being due to relative orientation of the disk and the line 
of sight - that BAL gas aff'ects an insignificant part of the 
outcoming radiation. 

As mentioned above, iV(Heii) of ~10^^ cm^^ corre- 
sponds to iV(Hi) of few times 10^^ cm~^ . This amount of 
neutral hydrogen gives origin to the continuum depression 
by ~15-30% starting somewhere above 912 A in the rest- 
frame - the shielding gas can not only outflow but infall 
as well (Murray et al. 1995; Proga et al. 2000) shifting the 
depression to the red. For non-BAL QSOs, the presence of 
such depression would indicate that there is enough Heii 
to produce a noticeable intensity break at energies above 4 
Ryd in the outcoming radiation. It seems natural to assume 
that the depth of this break depends on the angle to the ro- 
tational axis of the disk. A statistical study of the hydrogen 
Lyman continuum depression in the flux-calibrated high- 
resolution quasar spectra is expected to deliver a compre- 
hensive information about the covering factor of the Hell 



km s ^ ) absorption lines displaced by ~ 14000 km s ^ blueward 
of the quasar systemic velocity was detected at z = 3.021 to- 
wards CTQ 325 (Levshakov et al. 2004). 



opaque gas. This study is still to be done. However, there 
are several hints that this covering factor may indeed be 
signiflcant. Continuum depression starting roughly at ~930 
A (rest-frame) is clearly seen in HST spectra of non-BAL 
quasars shown in Binette et al. (2005) who state that such 
quasars (dubbed by them as 'class A' quasars) comprise 
about 60% of their sample (106 species). The estimates of 
the corresponding column densities of H l arc not very cer- 
tain due to the low quality of spectra, but nevertheless they 
vary between 2 and 6 x 10^^ cm~^ - well within the pre- 
dicted boundaries. A good example represents a HST spec- 
trum of HS 1103+6416 shown in Fig. 1 in Kohler et al. 
(1999) - continuum depression starting at ~950 A allows 
an accurate estimate of N{ILi) ~ 2.5 x 10^^ cm~^ . 

It is to note that so-called 'composite' spectra obtained 
by co-adding of many individual QSO spectra and shown 
in Fig. 5 in Telfer et al. (2002) and in Fig. 3 in Scott 
et al. (2004) are of no use in the present context since 
fine features present in an individual spectrum vanish or 
are smeared out in the co-added one. Moreover, before co- 
adding the individual spectra were 'statistically corrected' 
for the Ly-a forest absorption using the general distribu- 
tion of the Ly-Q clouds over A''(Hl) and z. This operation is 
quite subjective because statistical properties of the Ly-a 
clouds distribution vary significantly from sightline to sight- 
line. However, even in the composite spectra some feature 
are present which may be attributed to the quasar wind: 
these are a trough-like structure between 700 and 900 A in 
Fig. 5 (Telfer et al. 2002), and absorption features blueward 
of O VI and Ne viii emission lines again in Fig. 5 (Telfer et 
al.) and in Fig. 3 in Scott et al. (2004). The former may 
be caused by the Hi Lyman continuum absorption, and 
the latter - by the resonance absorptions in the Ovi and 
Neviii lines. 

An additional support for possible Heii opacity of 
quasar winds comes from X-ray observations which detect 
in ~50% of quasars so-called 'warm absorbers' - broad 

absorption feature due to the blend of various ionization 
edges, among them O Vii and O Viii (Piconcelh et al. 2005). 
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Fig. 12. Broad and shallow absorptions detected in the spectrum of HE 2347-4243. The zero radial velocity is fixed a,t z = 2.902. 
For other details see capture to Fig. 8. 



There arc evidences that these absorbers have A^(H) ~ 
10^2 — io23 cm"^ and the ionization parameters U < 5 
(Krongold et al. 2003; Schartel et al. 2005; Steenbrugge et 
al. 2005; Jinienez-Bailon et al. 2007). This leads again to 
iV(Heii) -10^* cm-2 . 

To summarize, it seems quite likely that the outcoming 
QSO radiation in general is not a pure power law but has 
an intensity break at £^ > 4 Ryd due to the Hell Lyman 
continuum absorption in the quasar disk wind. This can 
affect the rate of the He ii reionization in the intcrgalactic 
clouds. Recent observations of the Hi/He ii Lya forest to- 
wards HE 2347-4342 (Shull et al. 2004; Zheng et al. 2004) 
and HS 1700+6416 (Fechner et al. 2006; Fechner & Reimers 
2007) revealed significant fluctuations of rj which suggest a 
variable softness parameter S = JiR,yd/>/4Ryd of the meta- 
galactic radiation field. In this respect the revealed Hell 
opacity of the quasar wind which changes depending on 
angle to the rotational axis of the disk and by no doubt 
differs also from quasar to quasar can be responsible for 
the local fluctuations of S. For a pure power law ionizing 
spectrum, oc the photoionization rate of Heii is 
Fi cx 1/(q! + 3). For a spectrum with a break at 4 Ryd, 
tJ^°" > 0, the rate, becomes slower : 



, exp(-rH°"/x3) , , E 

r2 cx / TT^; dx. With X 



ra+4 



4 Ryd 



The ratio ri/F2 as a function of r^"" is shown in Fig. 14. 



3.2. Element abundances in quasar associated systems 

The properties of the NAL systems themselves are rather a 
by-product of the present study. However, there are some 
interesting results obtained which are worth to be men- 
tioned. 

Metal content as well as the relative elemental abun- 
dances (N/C, a-element/C, etc.) measured from the ab- 
sorption systems are usually used to conclude about the en- 
richment mechanisms of the absorbing gas. Concerning the 
quasar-related absorption systems analyzed in the present 
work, the measured ratio N/C show that these mechanisms 
can be different for gas in the vicinity of the central en- 
gine and in the host galaxy. For instance, the systems at 
Zabs = 2.198 (J 2233 606), z'abs = 2.107 (HE 1341 1020) and 
^abs = 2.352 (Q 0329 385), all with systemic velocities of 
thousands of km s~^ and, thus, parts of the outflows from 
the central region, reveal oversolar metallicity along with 
the relative overabundance of nitrogen as compared to so- 
lar values, [N/C] > 0. On the other hand, the systems at 
Zabs= 2.147 (HE 1341-1020) and Zabs = 2.898 (HE 2347- 
4243), with the systemic velocities of ~100-200 km s~^ and 
belonging probably to the gas in the host galaxy (cf. Nagao 
et al. 2006a), demonstrate similar oversolar metal content, 
but underabundant nitrogen, [N/C] < 0. Overabundance of 
nitrogen, a common finding also in the analysis of quasar 
broad emission Hues (e.g., Dietrich et al. 2003; Nagao et 
al. 2006b) , is supposed to be a consequence of rapid stellar 
evolution of quasar gas which leads to scaling of N/0 ratio 
with metallicity (Hamannfe Ferland, 1999). However, dif- 
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Fig. 13. Broad and shallow absorptions detected in the spectrum of Q 0329-385. The zero radial velocity is fixed at z 
For other details see capture to Fig. 8. 



2.435. 



ferent ratios [N/C] at similar mctallicity measured in the 
circumnuclear gas and in the host galaxy indicates different 
enrichment mechanisms which may be due to variations in 
the initial mass function. Note that underabundance of ni- 
trogen not complying with theoretical predictions of [N/C] 
and [Fe/C] was measured in the associated absorption sys- 
tems towards HE 0141-3932 (Reimers at al. 2005). 



4. Conclusions 

We have analyzed associated metal absorption systems 
identified in the UVES/VLT high resolution spectra 
{FWHM ~6 km s'^ ) of four quasars J 2233-606, HE 
1341-1020, HE 2347-4243, and Q 0329-385 with the aim to 
restore the spectral shape of the underlying ionizing contin- 
uum in the EUV range 1 < _E < 10 Ryd. The systems were 
selected on base of criteria which included the presence of 
lines of different ions and possibility to estimate accurate 
covering factors and column densities. This ensured the re- 
liability and robustness of the restored SEDs. All systems 
are physically related to the quasar/host galaxy and, thus. 



the SEDs in question are representative for the quasar out- 
coming radiation. 

The results obtained can be summarized as follows. 

1. The ionizing spectra responsible for the ionization struc- 
ture in the analyzed systems reveal a sharp intensity 
depression at _E > 4 Ryd which is attributed to the ab- 
sorption in the He II Lyman continuum. The required 
column density of Hen is of order 10^* cm~^ . 

2. The most probable source of the Heii opacity is the 
quasar accretion disk wind - a stratified fiow compris- 
ing dense near-equatorial gas accelerated by combined 
action of radiation-pressure and magneto-centrifugal 
forces (seen as UV BAL absorption) and more highly 
ionized shielding gas which protect the inner wind re- 
gions from ovcrionization by X-rays from the central 
source. All objects considered here do not belong to 
BAL-quasars which means that the He II opacity is due 
to the shielding gas. Assuming a pure power law SED 
at the wind base, we can estimate the column density 
of neutral hydrogen corresponding to A^(Hell) ^10"'^* 
cm~^ : iV(Hi) a few times 10^^ cm~^ . This amount 
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Fig. 14. The ratio of the photoionization rates Fi (a pure power 
law spectrum F^, oc to r2 (a power law spectrum with a 

break at 4 Ryd) as a function of r^**". 

of neutral hydrogen can be seen in a quasar spectrum 
as a weak continuum depression starting at A912 A 

(rcst-framc) and sometimes also as broad (stretching 
over 1000s km s~^ ) and shallow absorption in Hi Ly- 
a. In case of high metallicity (solar to oversolar) of the 
wind, the absorption features due to resonance lines of 
O VI in the FUV, Nc vi-Nc viii in the EUV, and O vii- 
O VIII in the soft X-ray ranges arc expected. In the UV 
spectra of quasars studied in the present paper, we do 
find broad and shallow Hi Ly-a and Ovi AA1031, 1037 
absorption features (Nevi-Neviii beyond the available 
wavelength coverage) as well as continuum depression 
at A ^912 A in the flux-calibrated spectra of J 2233- 
606 and HE 2347-4243 (flux-calibrated spectra of other 
analyzed quasars are not available). 

3. In order to estimate the fraction of quasar outcoming 
radiation affected by the Hen Lyman continuum ab- 
sorption, a systematic study of the predicted contin- 
uum depression in the flux-calibrated high-resolution 
QSO spectra at wavelengths shorter than 912 A (rest- 
frame) should be performed. This study is still to be 
done. However, observational data available now both 
in UV and X-ray regions suggest that at least ~50% of 
the quasar radiation passes through the gas opaque in 
the He ii Lyman continuum. This means, that the out- 
coming power law ionizing continuum has a pronounced 
intensity break at > 4 Ryd with a depth of this break 
depending on the angle to the rotational axis of the ac- 
cretion disk (the closer to the disk the deeper the break) . 
This can influence the rate of the He ii rcionization in 
the intergalactic medium and partly explain inhomoge- 
neous (patchy) ionization structure of the intergalactic 
Hell observed at z ~3. 

4. Due to careful selection and detailed analysis of the ab- 
sorption systems we obtained quite accurate estimates 



of the elemental abundances in the absorbing gas. It is 
shown that the relative abundance of nitrogen to car- 
bon, N/C, differs for the circumnuclear and interstellar 
gas of the host galaxy: at the similar oversolar metal- 
licity the circumnuclear gas is overabundant in nitro- 
gen, whereas the galactic gas is underabundant. This 
supposes different metal enrichment mechanisms which 
may be due to difl[erent IMF in the regions near the 
center and at the periphery of a quasar host galaxy. 
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